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Abstract 



We present some results about the interplay between the chiral and decon- 
finement phase transitions in parity-conserving QED^ (with flavours of 
massless 4 component fermions) at finite temperature. Following Grignani et 
al. 1 11, 12 1, confinement is discussed in terms of an effective Sine-Gordon the- 



ory for the timelike component of the gauge field Aq. But whereas in [11,12| 
the fermion mass m is a Lagrangian parameter, we consider the m = case 
and ask whether an effective S-G theory can again be derived with m replaced 
by the dynamically generated mass S which appears below Tc/j, the critical 
temperature for the chiral phase transition. The fermion and gauge sectors 
are strongly interdependent, but as a first approximation we decouple them 
by taking S to be a constant, depending only on the constant part Aq of 
the gauge field. We argue that the existence of a low-temperature confining 
phase may be associated with the generation of S; and that, analogously, the 
vanishing of S for T > T^h drives the system to its deconfining phase. The 
effect of the gauge field dynamics on mass generation is also indicated. 



*0n leave from Centro Atomico Bariloche, 8400 Bariloche, Argentina 
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I. INTRODUCTION 



Quantum electrodynamics in (2+l)-dimensions (QED3) is an interesting theoretical lab- 
oratory in which to explore, in a simpler environment, some of the fundamental features of 
more complicated four- dimensional gauge theories, such as chiral symmetry breaking and 
confinement; it may also be of direct physical relevance in some condensed matter systems 
|l|]. In this paper, we shall be concerned with A^- flavour QED3 at finite temperature, in 
its parity invariant (four-component) form. Our aim is two-fold: to cast some light on the 
nature of the phase transition associated with the dynamical generation of fermion mass 
(the "chiral phase transition"); and to argue for a dynamical interplay between the chiral 
phase transition and the confinement-deconfinement one. 

Chiral symmetry breaking ("csb"), or equivalently dynamical fermion mass generation, 
in QED3 at T 7^ has been extensively studied using Schwinger-Dyson equations, in the 
large N limit Though the approximations are admittedly not fully controlled, all 

calculations show that a fermion mass is dynamically generated for T < T^h, thereby spon- 
taneously breaking the chiral symmetry. But the Coleman-Mermin- Wagner theorem [|^,|^ 
forbids the occurrence of an order parameter transforming non-trivially under chiral trans- 
formations in this two-dimensional system. One presumes, therefore, that the transition at 
T = Tch is of the BKT type P,p!0| . However, to our knowledge there has been no explicit 
demonstration of this - for example, by mapping (or otherwise relating) the theory to one 
with a known BKT transition, which could then be identified with the csb transition in 
QED3. 

On the other hand, in a quite independent development, Grignani et al. [|n|'0 have 



shown that, when chiral symmetry is explicitly broken by a fermion mass m in the La- 
grangian, the relevant effective theory (at least in a certain region of parameter space) 
has a BKT transition at a temperature Tc, which the authors interpret as a confinement- 
deconfinement phase transition. 



It is natural to ask, first, whether the analysis of [pT| , p!2| can be extended to the case 
in which the fermion mass is dynamically generated, rather than appearing explicitly in the 
Lagrangian; and secondly, whether the resulting effective theory leads us any nearer to an 
understanding of the chiral phase transition, and of any possible connection between the 
two transitions. The purpose of this paper is to make a first attempt at answering these 
questions. 



It will be helpful to recall the argument of Grignani et al. |11,12|. According to them, 
confinement in QED3 at T 7^ can be probed by an Abelian analogue of the Polyakov loop 
operator, namely 

Pg(f) = e'^'/o''^o(^'^)'^^, (1) 

where (3 = 1/T [kB = 1), and e 7^ integer xe^ where ej is the charge of the dynamical 
fermions in the theory. The Euclidean path integral for the partition function is invariant 
under gauge transformations A^ = Afj_ + d^X: where 

X{(3,x) = x(0,f) + 27m/ Cf , n e Z. (2) 

The group of all gauge transformations modulo those which are strictly periodic is Z, the 
additive group of integers; this Z-symmetry is a global symmetry. Under a Z-transformation, 
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Pe(x) ^ Pe(x) e^™'^/"^ (3) 

If Z-symmetry is unbroken, then < Pe{x) >= and the system confines e charges; if 
Z-symmetry is spontaneously broken (as it can be, with an order parameter, being discrete), 
then < Pe{x) >7^ and the e charges are deconfined. 

As noted in [^,0, it is possible to choose a gauge in which Aq is independent of r. In 
this gauge, 

P~{x) = e^'^^^o^^) (4) 

and ^-transformations are 

efPAo{x) efPAo{x) + 27m, (5) 

under which Pg transforms as in @. The effective field theory for Pe{x) is then the effective 
action for the static field Aq{x), which is obtained by integrating out the other degrees of 
freedom from the path integral. The resulting effective action is in general non-local and 
non-polynomial, but it must still obey the Z-symmetry periodicity: 

S,ff[ef/3Aoix)] = S,ff[ef(3Aoix) + 2nn]. (6) 

A local approximation can be obtained by expanding Se/f in powers of derivatives of 
Aq{x) divided by the (Langrangian) fermion mass m. In such an expansion, the effective 
potential is obtained from the fermion determinant in the presence of a constant background 
field Aq - and it has to be periodic as in It is, one might think, very likely that such a 
potential should be essentially of "Sine-Gordon" type, and indeed this is what the authors 
of [0,0 were able to show, provided that 

m>r,ej/127r (7) 

(see also section II. A below). In this regime, the critical behaviour of the effective theory for 
Ao{x) is that of the 2-D Sine-Gordon model, which has a line of critical points along which 
a BKT transition occurs, separating confining and deconfining regions of parameter space. 
The critical temperature is 

e^/Svr 
~ l + e2/127rm + --- 

up to one-loop order. Thus the confinement-deconfinement transition has been shown to 
be of BKT type, by demonstrating that the relevant effective field theory has this critical 
behaviour. 

Consider now how dynamical (fermion) mass generation (dmg) might fit into this picture. 
Our basic idea is to see if it might be consistent for the dynamically generated mass E to 
play the role of the Lagrangian mass m in the analysis of [^,0. To investigate dmg, 
we shall use Schwinger-Dyson (SD) equations as previously, but here with one important 
difference: we wish to safeguard the crucial invariance (|^). We shall therefore split Aq into 
its slowly varying (Aq) and fluctuating (oq) parts, and - as a flrst approximation - solve the 
SD equations for the fermion and ao-fleld propagators in the presence of the uniform Aq 
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which, at this stage, plays the role of an external parameter, and which may be "large". 
The fermion self-energy equation will produce a dynamically generated mass S for every 
set of external parameters (in an appropriate region of parameter space) and for each Aq. 



The low momentum components of Aq are then treated as in |lTT| , pr2| with m replaced by E; 
for the analysis to go through as before we shall need to find that S does not depend too 
strongly on Aq, at least in some parameter region. This picture will only be consistent if 
(c.f. (^ ) the generated S satisfies the condition 

S>T,eJ/127r (9) 

with, of course, T < Tch- In this case, we can tentatively associate the generation of fermion 
mass with the occurrence of a confining phase. 

Of course, the fermion and gauge dynamics of Oq are coupled: the vacuum polarisation 
function will be the bubble diagram evaluated in the Ag-background, and using the dynami- 
cally generated mass E, which also depends on Aq. We shall see that in the low-temperature 
confined phase, where Aq is fluctuating about the trivial minimum of the Sine-Gordon po- 
tential, the properties of the gauge field propagator differ little from those in the normal 
{Aq = 0) case. However, for Aq configurations fluctuating near a position of unstable equi- 
librium for the S-G potential - which may be associated with the onset of the deconfining 
phase transition - we find that an instability develops in the vacuum polarisation, which 
threatens to destroy the non-zero gap solution and so restore chiral symmetry. There is an 
indication here, therefore, for an intimate connection between the chiral and the confining 
transitions. Unfortunately, our method does not allow us to explore the behaviour close to 
the phase transitions in any detail. 

The paper is organised as follows. In Section II we define the model and discuss the 
formalism used to study confinement and csb simultaneously. In Section III we calculate 
the vacuum polarisation in the presence of Aq. In Section IV we write down the self- 
consistent gap equation in the presence of Aq, and present numerical results together with 
our conclusions. An Appendix discusses the error due to neglecting the dynamics of the 
space components of the gauge field. 

II. TOWARDS A UNIFIED TREATMENT OF THE CHIRAL AND 
DECONFINEMENT PHASE TRANSITIONS 

In this work we shall be concerned with a theory whose generating functional of complete 
Green's functions fj, rj], in the imaginary time (Matsubara) formalism, is defined by 



Z[j^-f],rf^ = J VA^V^V^ exp^-S + J d'x[j^{x)A^{x)+r]{x)^{x) + tlj{x)7]{x)]\ (10) 
where / d'^x ■ ■ ■ = Jq dr J d'^x ■ ■ ■, /3 = 4, and S denotes the Euclidean action 



S 



j d^xC C = Cg + Cf ■ (11) 



The gauge field Euclidean Lagrangian, Cq, is assumed to be of the standard Maxwell type. 
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= \ F^^uF^u F^, = d^A, - d,A^ , (12) 

where we have omitted, for the time being, the gauge fixing. The fermionic field ip denotes 
fiavours ipa-, each one minimally coupled to the gauge field A^. The corresponding Euclidean 
Lagrangian, Cp, is then 

TV 
a=l 

where the coupling constant = (replacing the e/ of Section I) has the proper A^- 
dependence to ensure that, in the large- limit keeping e fixed, all Feynman diagrams are 
finite. The Lagrangian mass m may or may not be present. We use the notation r and 
X to refer to the time and space coordinates, respectively. The 7-matrices are in a 4 x 4 
(reducible) representation of the Dirac algebra in 2 + 1 dimensions: 




7. = I I ' (14) 

where /x = 0, 1, 2, and 

^0 = (I °i ) on = ( T n ) ^2 = ( " n" 1 ' (15) 



















(T2 = 





are the usual Pauli matrices. This representation allows for the introduction of a parity- 
conserving mass term for the four-component fermions. This corresponds, in the two- 
component formalism, to mass terms with opposite signs for each fermion field within a 
given fiavour. We shall here use a description that allows us to study the possibility of both 
chiral and deconfinement phase transitions simultaneously. Note that the latter transition, 
as presented in [ pT] , pr2[ , is considered in terms of an effective action for Aq, the time com- 
ponent of the gauge field. This effective action is derived in the static-Ao gauge, performing 
a derivative expansion. The usual way to study the chiral phase transition is, on the other 
hand, to work with truncated Schwinger- Dyson equations (SDE's), an approximation which 
is justified in the large-iV limit. It is possible, as we shall see, to find a non-empty region 
in parameter space where the simplifying assumptions made for the study of each phase 
transition are indeed compatible. 

We shall briefiy review here the derivation of the effective action for the field Aq, for the 



case of massive fermions, following |]TT],|I2[. This effective action is the interesting object 
to study when dealing with the deconfinement phase transition. We will afterwards extend 
this procedure, in order to include the non-perturbative dynamics responsible for the chiral 
phase transition. The fermion mass, rather than being a Lagrangian parameter, will then 



be dynamically generated, and determined by a gap equation. Following references p!l|JT2 
(with trivial changes due to the fact that we now have N fermionic fiavours) the effective 
action S'e//[Ao] for a system where the fermions do have a Lagrangian mass m, is obtained 
by first using a gauge transformation such that Aq is r-independent [|, and then integrating 
out all the remaining fields, i.e.. 



^Aq itself cannot be set equal to zero, because the time coordinate is compact. 
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where the functional rAr[v4o, Aj] is the resuh of integrating out the fermionic fields, 

= {det[^ + ieN'yoAQ{x) + ieN'yjAj{T,x) + m]}^ . (17) 

Of course, we have rAr[y4o,Aj] = A^r[ylo,Aj], where r[y4o,Aj] denotes the contribution 
corresponding to the integration of just one fermionic flavour. Following [11,12|, r[y4o, Aj] 
may be evaluated by using a derivative expansion for the external gauge field, taking into 
account the fact that Aq may have 'large' values, i.e., its constant component cannot be 
assumed to be small, since that would spuriously break the invariance under shifts in Aq. 
To use an approximation which preserves this periodicity is crucial for the study of the 
deconfinement phase transition. Thus we let 

A^{x) = A^{x)+a^{x) (18) 

where Aq is the large, almost constant piece, Aj = 0, and a^{x) denotes the fluctuating part 
of Af^{x). As Aq cannot be assumed to be small, the expansion of the determinant must 
include that object exactly: 

det[^ + ieN'yoAo{x) + ieNljAj{T,x) + m] = det[^ + ie^^yoAo + m] 



X det 



1 + icNi^ + ieNloAo + m) ^ (^(t, 5) , (19) 



where, for the last factor in ([TO), we can attempt a small-a^ expansion. In a first approxi- 
mation, and to be consistent with the approximation usually invoked in the SDE's approach 
to the study of the chiral phase transition, we shall ignore the dynamics of aj, the space 
components of the gauge field. A discussion about the error due to this approximation is 
postponed to the Appendix. Then, we only need the fermionic determinant as a function of 
the time component of the gauge field. For one flavour, we have 

g-r[Ao] ^ det[^ + icN-foM^) + m] 
= det[^ + icN'yo^o + m] det 1 + ieN{^ + ie^'joAo + my^'yoao{x) . (20) 



The ao-independent factor in (|T9D, involving no derivatives of the gauge field, plays the role 
of an effective potential, and is given by 

det[^ + ieAr7oio + m] = e/'^'-^(-.^) (21) 

with 

V{m,^) = ^logdet[{-zdo-er,Aof-V' + m']. (22) 
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Regarding the second factor, when expanded for low momentum, it produces a wave function 
renormahzation of the kinetic term for Aq. Putting together these two pieces, and including 
the factor A^, we see that the effective action for Aq becomes 

Ses[Ao] = J d^x (^Z(m,e^io/T)^VAo ■ VAq - iV V^(m, e^v^/T)) , (23) 

where 

ivnoo(o,P) = iVnoo(0,0) + P(Z-l) + ... (24) 

where noo(^, k"^) denotes the 00 component of the vacuum polarization tensor in the presence 
of the constant gauge field Aq. 

To relate this form of the effective action to a Sine-Gordon description, one performs a 
harmonic expansion of the effective potential: 

Vim,ej,A/T) = ti-ir e— (l + ^) . (25) 

This effective potential may be further simplified in some particular regimes. For exam- 
ple, when m » T, T/m and e^/(127rm) are small and e^/T is arbitrary, so that Z \ 
and the higher harmonics in the effective potential are exponentially small perturbations; 
one may then keep only the leading term, 

/ 7T7 

V{m,e^A^lT) e— cos(e;v^/T) , (26) 

vr 

which is a sine-Gordon potential. Under these conditions, one has an effective description 
which corresponds to the Sine-Gordon theory for the scalar field Aq : 

^ir^^^[^o] = / dx\^^\/A, . VAo - iV^ e-'"/^cos(e^Ao/T)} , (27) 

which has a BKT P,p!o| phase transition. The "/?" parameter of the Sine-Gordon model is 
now dependent, and given by /5 = Under the assumption (|^ (with e/ replaced by 

Cat), the critical temperature is given by (^) with e/ replaced by cat. 

It is important to remark that, in this analysis, the Lagrangian mass m has been regarded 
as a free parameter. We shall now argue that, when m is zero, one can - under certain 
conditions - arrive at an effective theory of the same form as (pTf), but with m replaced by 
S, the dynamically induced gap. 

The effective action (|27|) has been obtained by performing a derivative expansion. In 
consequence, when the Aq field is regarded as dynamical, the effects due to its large mo- 
mentum modes are either neglected, or improperly taken into account. This is, of course, 
not an issue when one is dealing with the deconfinement phase transition, but it will have 
an important influence for the particular case of massless fermions. In this case the mass 
parameter used to perform the low-momentum expansion is proportional to the tempera- 
ture, and one is then missing modes of the gauge field which may indeed be responsible 
for the generation of a dynamical mass for the fermions. In other words, the procedure 
of using the low-momentum effective action for Aq for the case of massless fermions would 
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not be accurate, since one still has to include virtual corrections which contribute to the 
generation of a mass for the fermions. If a non- vanishing dynamical mass for the fermion is 
generated, the momentum expansion procedure can be made consistent, but one should be 
careful to include in that dynamical mass the dependence on the temperature and on the 
parameters of the model. 

We shall now give a formal derivation to justify this procedure. The definition of the 
effective action for the large field Aq should now take into account the existence of the high 
momentum modes for A^. The natural way to proceed is to decompose the gauge field into 
its large and small components, as in (|TBp, where only up to two derivatives of A will be 
kept, while is assumed to take care of the higher derivatives. We shall, to simplify the 
notation, get rid of the tilde, and A will denote just the large component of the gauge field. 
Then, Se/f is defined by 

exp{-Seff[Ao]} = J T)a^Ila'D'ipa'Dipaeyip{-S[Ao + ao, aj,ijj,ip]} (28) 

where S is the Euclidean action, as defined in (|ll]). More explicitly. 



S[Aq + oq, aj, i/j^ip] = / d^x 



^F^,{A + a)F^,{A + a) + Y,M^ + ^eN{4+ 



■ (29) 



It is crucial to realise that A and a are assumed to have complementary supports in mo- 
mentum space. Thus the mixed quadratic term for the gauge field in (p9| ) vanishes 



Jd'x^F^,{A)F^,{a) = 0, (30) 



as can be easily verified by Fourier transformation. 
Then we see that 



exp{-Seff[Ao]} = exp{-i J d^x{d,Ao)' - W[Ao]} (31) 



where 



^-wiAo] ^ / Va^Vi,V^lJexp{- / d'x[-F^,{a)F^,{a) + + teNi4+ ■ (32) 

Note that Ty[v4o] has a simple interpretation: it is the free energy corresponding to parity 
conserving QED^ in the presence of an Aq background. The leading term in a large N 
expansion is 

W[Ao] = -TT[\nSp'] (33) 



where Sp denotes the full fermion propagator, for the theory defined by ([3^). This propa- 
gator is the result of solving the S-D equations in the Aq background. 

In principle, since Aq depends on x so will S, but we shall make the assumption that the 
induced mass is not sensitive to the gradient of Aq, but only to the constant part Aq. This 
assumption will only be reliable in the confined phase, where the topological excitations 
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of the S-G theory are absent, and Aq has only small fluctuations around one of the stable 
minima of the effective potential. With E assumed constant, Sp'^ will have the form 

Sp^ =^ + teNloAo + J:; (34) 

evaluation of the trace in ( ^31) may then proceed as in the m 7^ case, via Equations (^OD 



(P^) above. But an important question to answer will be whether, for at least some part 



of the solution-space for E, conditions (|^) do in fact hold, allowing the passage from ( pSf ) to 



The low-momentum effective action for the gauge field is then constructed with the IPI 
functions derived (see Sections III and IV) from the self-consistent SD equations, with Aq 
appearing as an external field. As the low-momentum component of the gauge field has 
not yet been integrated, one has to treat its dynamics separately; this is done with the 
S-G description, as in |TT],|T^, which does indeed allow for a BKT transition. The phase 



transition for the deconfinement transition should then be studied taking into account the 
fact that S is now not an independent parameter (as was m), but rather a function of T, e, 
N and Aq. In this connection, we shall have to check that the dependence of S on Aq is, 
in fact, not so rapid as to change the effective potential in ( [27| ) significantly away from the 
S-G form. 



III. THE VACUUM POLARISATION FUNCTION nnn IN THE PRESENCE OF i. 



We derive here the vacuum polarisation function in the one-loop (large- A^) approxima- 
tion, keeping the full dependence on the 'large' (uniform) piece of the time-like component 
of the gauge field, Aq. 

In principle, the coupled S-D equations for the IPI functions should be solved self- 



consistently, having made the usual truncations. However, experience shows [jT3[ that, 
in dynamical mass calculations, it makes very little difference whether a fermion mass is 
included in IIoo or not. We shall present the calculation for massless fermions, indicating at 
the end the changes needed to deal with massive ones. 

For one single four-component fermionic flavour (for flavours we only have to multiply 
the result by A^), the object we consider is the quadratic form r(^)[Ao], deflned by 

r('^[^o] = ^ / ^So(0,-A;)noo(0,A;)5o(0,A;) (35) 

where a^(n, k) denotes the Fourier transformed of a^, namely 

1 +00 

(^,ir,x) = -5 E e'^'^'^^^'-^^~a,{n,k) (36) 

P n=—oo 

with an = the bosonic Matsubara frequency. From the deflnition of the fermionic 
determinant, we see that the vacuum polarization tensor is given by 
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(=+oo 

E 



(2^ 



tr 



—Iv 



(37) 



where we adopted the notation uji = uji + e^Ao, where ui = {21 + 1)^ is the fermionic 
Matsubara frequency. Note that the Dirac slash is now defined for the two spacehke indices 
in the 'dimensionally reduced' theory. The 'tr' denotes Dirac trace, over the four component 
spinor space. After evaluating this Dirac trace, and performing the angular integration for 
the spatial momentum, F{k) = noo(0, k), the zero frequency part of IIoo may be written as 
follows 



-N 



da 



dp^ 1 



TT 



+ 00 

E 



p^ — a{l — a)k^ — Lul 



P ^^^[p' + ail-a)k^ + u, 



(38) 



where a is a Feynman parameter. To calculate F{k), we found it convenient to split that 
function into FT=o{k) (its T = limit) and F{k) = F{k) — FT=o{k). By performing the sum 
over the discrete frequencies and the integration over the modulus of the spatial momentum, 
the subtracted function F{k) may be expressed as an integral over a: 



F{k) = % /'c/aln[l + 2cos(ejvMo)e-^v/^(^'= + e-2'3y^(i^'=] 



(39) 



131 , and given by a^k/S. 



The zero-temperature value of F is well-known 

It does not seem possible to evaluate ( p^ ) in closed form. However, since IIoo appears 
inside an integral in the equation for the dynamically generated fermion mass (Eq.(^5D 
below), it would make the solution of ( p5D much easier if an accurate analytic approximation 
were available for (^). A similar problem was encountered in 0, and there a convenient 
approximation was found which preserved correctly the limits — > and (3k ^ oo (see 
Eq.(9) of PI). Following the same procedure here, we approximate IIoo by 



which incorporates 



Tif3k 



TT 



81n[2 + 2cos(ejv/3Ao)] 



(40) 



and 



Jim noo(0,A;) 

pfe— »oo 



e'^k 



(41) 



cos 



\ 2 , 



(42) 



and where we have now written the result explicitly for 7^ 1 fiavours (recall that 



e/yN, and e is held fixed as A^ varies). Expression ( ^2]) is the "thermal mass", in the 
presence of Aq. It is clear that when c^I^Aq > ^ the thermal mass goes negative. We 
shall discuss this further in Section IV. For the moment, note that this value is in any case 
approaching near to the "unstable vacuum" of the S-G potential at cn^Aq = vr, which lies 
beyond the region of validity of our approach, based as it is on perturbing around one of 
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the stable minima of the S-G potential. For cnPAq ~ 1.5, we find that the approximation 
(^) is always within 10% of the exact formula (0). 

It is not hard to show that, when m ^ 0, the zero frequency part values of IIoo are given 
by a function F{k) defined by 

Fik) = I da I - - [,2 + ,(i_,)fc2 + ^2 + ^2]2 ' (43) 

An analogous procedure to the one followed in the massless fermion case allows us to 
write the subtracted function F as 



2 f) 1 

F(k) = ^ (1 - m^) / da ln[l + 2 cos(e^/3io)e-^V^(^^^)^^ + e-2^V°(i-°)^''+™'l . 
7r/3 am Jo 

(44) 



As for the massless case, an approximate expression for F{k) may also be obtained, but 
in our calculations we have not used the massive form of IIoo. In a fully self-consistent 
treatment, m would be replaced by S in (^4]). 



IV. DYNAMICAL MASS GENERATION IN THE PRESENCE OF Aq 



In this first exploration of the interplay between csb and confinement in QED3 at T 7^ 0, 
we shall adopt the simplest approach to the S-D equation for E (the "gap equation" ) which 
was taken in the early paper by Dorey and Mavromatos . We shall assume that S is 
momentum-independent, that fermion wave function renormalization may be neglected, and 
that only the instantaneous-exchange part of the kernel is retained (see for evidence 

that relaxing these assumptions will not change the qualitative conclusions very much). The 
new feature of our calculation of S, as for IIoo, is the inclusion of the Aq background. 

The fermion propagator has the form (z + i'joLJn)'^ where u;„ = Un + cnAq and Un is 
the fermionic Matsubara frequency as before; the gap equation is then 



1 



kdk 



AiiNJ [F + noo(fc,AAo)]v/FTs^ 



sinh/3v/FTS2 



cosh /5 VFTS2 + cos(e/3Ao/ ViV) . 



(45) 



where IIoo is the /x = z/ = component of the vacuum polarization tensor 11^,^ as calculated 
in the preceding Section, Eq. (0), and is fixed as varies. Eq.(^5D is to be compared 
with Eq.(20) of 0, to which it reduces when Ao ^ (note that the latter equation is in 
error by a factor of | on the right hand side: "2A^7r" should read "4A^7r", so that the N = 2 
results of [0] correspond to our = 1 results, etc.; note also that our IIoo is equal to k"^ x IIoo 
of 01, but is the same as IIq of ||^). The upper limit of integration in (^5]) is taken to be 
unity, as in 0. 

It is worth noting that the momentum integral in the S-D equation for S should really 
have an IR cutoff e, given by a small fraction of S, since the gauge field corresponds to 
the fluctuating part of Aq. This should not make any difference to the solution of the 
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equation, for all the cases where there is a non- vanishing S, since there the IR behaviour of 
the integrand is smooth. We have verified numerically that indeed the solution to (^) is 
insensitive to the introduction of a small infrared cutoff. 

On the other hand, when the only stable solution is S = 0, the IR cutoff would make a 
difference, but then the separation into low and high momentum modes collapses. Again, this 
would correspond to the physics near to the transition, where our procedure is not applicable. 

We turn now to a discussion of the numerical solutions of (^), and their implications. 
Figure 1 shows s = ^ as a function oit = ^ for various values of 6 = — , for the case = 1. 
For 6 = 0, the curve reproduces the known results of 0. For 6 7^ 0, a mass continues to be 
generated; indeed, at a given temperature the mass increases with b (see Figure 2 where s 
is shown as a function of b for fixed t). However, for 6 7^ we cannot probe temperatures 
below 1^, since below this value the thermal mass goes negative, as remarked after (^2]). In 
fact, from ( ^0]) it is clear that when t = ^ (i.e. epAo = ^), IIoo reduces to the zero 
temperature value. This is the reason why the curves all tend to s{t = 0) as t tends to 

J. _ 3b 
27r • 

In reality, the calculation should probably only be accepted for t > - (i.e. cPAq < tt). 



since we must remember that, at cPAq = vr, the Sine-Gordon potential in ( pTl) has a point of 
unstable equilibrium, while our mass calculation is understood as being performed around 
a stable minimum of the potential. Inspection of Figure 1 shows that the condition t > ^ 
is met for the curves corresponding to 6 = 0.01 and 0.02, but only for part of the b = 0.03 
curve, and for less of the 6 = 0.04 curve. 

Figure 2 shows s as a function of b for various values of t, also for the case = 1. As 
before, the curves tend to the zero-temperature value as b tends to ^t. For the values of t 
and b shown, the more conservative limit b < nt is always satisfied. The critical temperature 
for 6 = is near tch = 0.0137, which is why the curve for this value of t "divides" the curves 
in Figure 2. 

We must now consider whether these results support our approach a posteriori. First 
and most importantly, we need to check whether conditions (^) hold. Inspection of Figure 1 
shows that indeed the condition E >> T (i.e. s » t) is met for small T, say T < 0.01. The 
condition s >> is less satisfactorily met, but it is not flagrantly violated. Secondly, we 
see from Figure 2 that, while there is certainly a dependence of S on Aq, it does not seem to 
be sufficiently rapid to call into question the effective S-G description of the deconfinement 
transition. 

So far we have discussed the case A^ = 1, and we must now remember that some ap- 
proximations we have made depend on the large- A^ limit. However, it is well known (see 
for example and references therein) that in the model considered here there is a critical 
value of A^, Ag(T), above which chiral symmetry is restored. Ag is typically in the region 
1-3. We shall therefore consider the case N = 2 for definiteness. While A^ = 2 is hardly a 
large value, there is evidence to suggest that the inclusion of terms which are higher order 



in 1/A^ will not alter the qualitative conclusions |14 . 

In Figure 3 we present, for N = 2 , the curves of s as a function of t, for two values of 
b. We see that both the critical temperature and the value of s at t = are reduced, but 
still there is an important fraction of parameter space where s » t. As for the condition 
s >> ^2nN ' satisfied to the same extent as for the A^ = 1 case, since going from A^ = 1 
to A^ = 2 reduces both sides of that inequality by (roughly) a factor 2. 
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We may conclude that we have demonstrated the feasibihty of our approach, in which 
the dynamically generated S replaces (at least in a limited region of parameter space) the 
Lagrangian mass m of [jTT],|12[. We may therefore argue that in the m = case, dynamical 



generation of S is associated with a low-temperature confining phase, via the emergence of 
an effective S-G theory. There is also an indication of a (reciprocal) effect of the gauge field 
dynamics on mass generation. As noted after Eq. (|42D , when cn^Aq > (27r/3) the thermal 
mass goes negative, which would imply a singularity in (^) due to the ensuing pole in the 
gauge field propagator in the Euclidean region. We may interpret this as indicating that, as 
the gauge field fluctuates towards the value cnPAq = n, which is an unstable stationary point 
of the S-G potential, the stability of the gap equation is lost, and mass generation ceases. 
Interestingly, apparently similar special temperatures occur in a study of the Gross-Neveu 



model with an imaginary chemical potential |T5[. Clearly it would be desirable to consider 



more fully the coupled dynamics of the Aq and fermion sectors, including an x-dependence 
for E. Nevertheless, as a first orientation, our procedure seems to be reasonably justified. 
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APPENDIX : INCLUSION OF THE SPATIAL COMPONENTS OF THE GAUGE 

FIELD. 



In the SDE's approach, one usually neglects the dynamics of the spatial components of 
the gauge field. This produces an error in the result of the self-consistent gap equation. In 
the evaluation of the effective action for Aq, on the other hand, the spatial components of the 
vacuum polarization tensor yield a wave function renormalization for the spatial components 



of A^. 



Let us call SSefj the contribution arising from the integration over the spatial modes a 
in (]TB|). It is quite straightforward to see that 



JVaex.pl- J2 iJ J^Pjki-n,-k)Fjk{n,k) 



( +00 \ r cPk 1 

X expj- ^ -y ^^--y^ai(-n, -A;)n,j(n, fc)aj(n, (46) 

where the tilde denotes Fourier transform. In the static Aq gauge there is a remaining gauge 
invariance under space dependent gauge transformations, which implies that the function 
H-ij is spatially transverse: 

Ilij{n,k) = U{n,k){-A6ij + didj) (47) 

Fixing these gauge field components to be spatially transverse, 

r { 1 r d'^k ~ 1 

e-'^^ff = JVaexpl- ^ - j -^^F,,{-n, -k)[l + U{n,k)]F,,{n,k)\ . (48) 

Then we may write for the correction to the effective action 

^-ss^ff = n„.[i + n(n,fc)]-i (49) 



or, equivalently, that there is a correction to the effective potential V, namely, V V + 6V, 
with 

5y = -X:/(01n[l + n(n,^)]. (50) 

In spite of the fact that this contribution exists, it is negligible in the large limit (note 
the absence of the factor which appears when integrating fermions). Moreover, even for 
only one fermionic flavour, we checked that 6V is exponentially suppressed in comparison 
with V. 
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Figure 3 
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FIGURE CAPTIONS 



• Figure 1: Curves of s = ^ as a function of t 
the case = 1. 

• Figure 2: Curves of s = ^ as a function of b 
the case N = 1. 

• Figure 3: Curves of s = ^ as a function of t 
case N = 2. 



for various values of 6 = — , for 
— for various values of i = ^, for 
^ for two values of 6 = for the 
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